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Abstract We combine Halo/Cluster Effective Field Theory (H/CEFT) and the Gamow Shell Model 
(GSM) to describe the + ground state of 6 He as a three-body halo system. We use two-body interac- 
tions for the neutron-alpha particle and two-neutron pairs obtained from H/CEFT at leading order, 
with parameters determined from scattering in the P3/2 and So channels, respectively. The three-body 
dynamics of the system is solved using the GSM formalism, where the continuum states are incorpo- 
rated in the shell model valence space. We find that in the absence of three-body forces the system 
collapses, since the binding energy of the ground state diverges as cutoffs are increased. We show 
that addition at leading order of a three-body force with a single parameter is sufficient for proper 
renormalization and to fix the binding energy to its experimental value. 

Keywords Nuclear effective field theory • Gamow shell model • Halo nuclei 



1 Introduction 

Nuclei located far away from the valley of /3-stability display peculiar features that do not occur for 
well bound nuclei. The strong coupling to the continuum manifests itself in the existence of halo 
configurations, where some nucleons orbit far away from a core of more tightly bound nucleons, and 
of Borromean systems, where removal of one nucleon is accompanied by at least one more nucleon. 
The neutron-rich Helium isotopes 6 He and 8 He offer two examples of such nuclei: both are Borromean 
halos that have no bound excited states. They also exhibit the "binding-energy anomaly", i.e., higher 
one- and two-neutron emission thresholds in 8 He than in 6 He. 

Halo configurations are characterized by a large nuclear radius compared to the size of the tightly 
bound core or, equivalently, by a small nucleon separation energy compared to the core binding energy. 
The physics of halo nuclei is a perfect arena for the application of effective field theory (EFT). EFTs 
provide a powerful framework to exploit separation of scales in physical systems in order to perform 
systematic, model-independent calculations. If, for example, the relative momentum k of two particles 
is much smaller than the inverse range of their interaction, l/R, using contact interactions observables 
can be expanded in powers of kR [lj. 

The application of EFT to halo and cluster systems, Halo/Cluster EFT (H/CEFT), was first 
exemplified in low-energy neutron-alpha particle (na) scattering 0. Even though there is no bound 
state, the na T matrix has a resonance pole at an energy E§ gs ~ 0.8 MeV much smaller than the 
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E^ex — 20 MeV it takes to excite the alpha particle. The physics in the p wave is closely related to halo 
dynamics, because the ground state in 5 He can in a first approximation be described as an na system 
in the p 3 / 2 -wave configuration, which has the qualitative characteristics of a two-body halo nucleus. 
H/CEFT captures these features at leading order (LO), and provides at next-to-leading order a good 
description of na scattering H/CEFT has also been successfully applied to other dilute two-body 
systems such as low-energy aa Q and proton- 7 Li Q scattering, radiative neutron capture on 7 Li @, 
and the electric properties of 11 Be H. 

The natural follow up to Ref. [2| is to consider the next halo system within the Helium isotope 
family, that is, 6 He. This is the aim of the present paper. Since the 6 He ground state is bound by 
only Eq 9S ~ — 1.0 MeV, it is appropriate to assume that it can be described as the three-body system 
n + n + a. The na interaction is that studied in Ref. |2j, while the nn force is a contact interaction 
determined by the low-energy nn scattering parameters 0; • 

It is well understood that the physics of three-body systems can be much richer than the physics 
of its two-body subsystems. A famous example is the Efimov effect which occurs for non-relativistic 
particles with short-range interactions: if the s-wave scattering lengths of its subsystems are tuned 
to infinity, there can be an infinite sequence of three-body bound states that has an accumulation 
point at the three-body threshold |8j. A closely related phenomenon is the Thomas effect [9j, where 
a finite-range two-body potential that is only attractive enough to support a single two-body bound 
state can produce three-body bound states with arbitrarily large binding energies as the range goes 
to zero. In EFT at a given order, the Thomas effect is a consequence of an inappro pria te omission of 
a three-body force. As it has been demonstrated for s-wave two-body interactions [1(|, a three-body 
force is necessary and sufficient for renormalization-group (RG) invariance at LO, allowing three-body 
energies to be independent of the ultraviolet regulator. The parameter associated with this force then 
provides a scale for the remaining discrete scale invariance, which reflects itself in the Efimov spectrum. 
The structure of 22 C and other two-neutron halos with s-wave neutron-core interactions was discussed 
using H/CEFT in Ref. [u|. 

The situation with p-wave interactions is less clear. There is debate over whether the Efimov effect 
can be realized in this case — see Ref. [12J and references therein. The main issue to be addressed 
below is whether a three-body force is needed at LO so that the EFT description of 6 He is properly 
renormalized. 

A microscopic description of weakly bound/unbound nuclei requires taking into account the inter- 
play between bound states, scattering states, and resonances. In other words, these systems have to 
be described as open quantum systems (OQSs), in contradistinction with well-bound nuclei, which are 
nearly isolated from the environment of scattering states and decay channels ("closed quantum sys- 
tems"). A recent realization of the shell model for OQSs is the so-called Gamow Shell Model (GSM) 

El m EI M- 

The GSM is based on the Berggren basis [17| . which consists of bound, resonant 
and scattering single-particle wave functions generated by a finite-depth potential, and provides the 
mathematical foundation for unifying bound and resonant states — the poles of the T matrix — in the 
context of the Schrodinger equation. The GSM has been used [1J; [l5|; llfj to study the properties of the 
Helium isotope family using a phenomenological na potential and a residual two-neutron interaction 
adjusted to few-body energies. 

In this paper we use the formalism of the GSM to solve the Schrodinger equation describing the 
dynamics of 6 He with the contact interactions of H/CEFT, including a possible three-body force. We 
focus on the binding energy of the + ground state. 

The paper is organized as follows. In Sec. [5] we review the potentials derived with EFT and used 
for the study of 6 He. We introduce the GSM formalism in Sec. [3J and in Sec. 2] we show results for the 
ground-state energy of 6 He. We shall see that without a three-body force, the system is not properly 
renormalized at LO Q- We conclude and summarize in Sec. [5l 



2 Two-Body Potentials 

Being weakly bound compared to the first excited state of the alpha particle, the — + ground state 
in 6 He can be described as a three-body system n + n + a. The neutrons in the halo interact with the 



1 Our first resu lts were presented in Ref. [T||- Similar results have been obtained independently by Ji, Elster, 
and Phillips [13 . 
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alpha particle via a two-body interaction V na and among themselves via a potential V nn . We denote 
the neutron (core) mass by m n (M c ) and the neutron-core reduced mass by fi = m n M c / (m n + M c ). 

The potential between the a core and a neutron is constructed with EFT as described in Ref. 

The small relative momentum means that neutron and alpha particle see each other, in a first 
approximation, as elementary objects. At LO there is only one contribution, which is in the p 3 / 2 
channel, and the "dimeron" potential projected onto this channel can be written as 



V na (k\k,k ) = 



k'k 



A 



Bkl 



(1) 



where k (k') is the incoming (outgoing) relative momentum and fco = \J1[iF na in terms of the to- 
tal energy E na of the na subsystem. A and B are parameters. Since this interaction is singular, a 
regularization procedure is introduced in form of an ultraviolet cutoff A na . The cutoff separates the 
short-distance physics, which is not included explicitly in the dynamics at low energies, and the long- 
distance physics, which is. This is here achieved by introducing a smooth regulator function 



F(x) — exp(— x), 



(2) 



whose role is to suppress the high-energy contributions of the potential. We thus replace the potential 
©by 

V na (k', k, fc ; A na ) = — A JOL- ^ F (fc' 2 ML) F (fc 2 ML) ■ (3) 



A(A r , 



B(A na )k 2 



In order for observables to be RG invariant, i.e., independent of the arbitrary cutoff, the parameters 
A(A na ) and B(A na ) must depend on A na . 

More precisely, at LO, A(A na ) and B(A na ) are fixed such that the phase shifts at low energies 
obtained with the potential ((3]) reproduce the effective range expansion (ERE) in the P3/2 channel 
truncated at the level of the effective "range" : 



fc 3 coU na (fc) = -— + ^ 2 , 



(4) 



with the scattering volume a na = —62.951 fm 3 and the effective momentum r na 
The position k res of the p 3 / 2 resonance is obtained from 

COt &na {kres ) ^, 



-0.8819 far 



(5) 



so at this order k res = (0.174824 — 0.031319i) fm 1 . By solving the Lippman-Schwinger equation with 
the potential ([3j), one obtains 



A(A na ) = 2 M 
B{A na ) = -n 



1 



A :i 



4V27T 

4 



2 

7T 



(6) 
(7) 



The two neutrons in the halo have sufficiently low relative momentum that meson exchange can 
be considered a short-range force. The neutron-neutron potential is thus taken from the pionless EFT 
At LO, the potential is entirely in the 1 So channel; in momentum space it is simply a constant 
C. As before, the potential requires regularization, for which we continue to use the function F(x), 
but now in terms of the relative momentum between the two neutrons and of an nn cutoff A nn : 



As previously, we fix the coupling constant C(A 
at the level of the scattering length, 



V nn (k', fc; A nn ) = C(A nn ) F (k' 2 /A 2 nn ) F (k 2 /A 2 nn ) . (8) 

with the ERE for nn scattering, but now truncated 



fccot^ nn (fc) = -- 



(9) 
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with a nn = —18.7 fm 21]. Again solving the Lippman-Schwinger equation, 



C(A nn ) = — 



1 



Ar. 



(10) 



Note that we do not modify the nn potential in 6 He to account for the presence of the a core, as 
frequently done [2 [l5| • This modification is a three-body effect that in EFT is represented by three- 
body forces, which are present starting at some order, since they are not forbidden by any symmetry. 
We want to determine whether such a force is needed at LO to renormalize the n + n + a system. 



3 Schrodinger Equation with the Gamow Shell Model 

We now consider the solution of the Schrodinger equation for the n + n + a system with the Gamow 
Shell Model. We use coordinates inspired by the Cluster Orbital Shell Model [12; HH : r $ is the position 
of neutron i = 1,2 relative to the a core, and Pi the corresponding momentum. The Hamiltonian of 
the n + n + a system with the two-body interactions V na and V nn is written as 



H 



E 



+ V nn (A nn ) + 



Pi • P2 

Mr ' 



(11) 



This Hamiltonian is translationally invariant, the recoil term pi • p%jM c stemming from the choice of 
coordinates. 

We work within the framework of the Gamow Shell Model formalism Q OH (HI to solve the 
dynamics generated by the Hamiltonian (|11[) . The three-body equation is solved using a single-particle 
(sp) basis. The set of sp states that define the one-body valence space is taken as the set of eigenstates 
of the LO potential V na (ko; A na ). They are solutions of the one-body Schrodinger equation 



H sp \W) 



Fp2 

1\i 



+ V na (k ; Ana) 



By inserting the completeness relation projected on the p 3 / 2 partial wave, 



dkk 2 \k)(k\ = 1, 



(12) 



(13) 



along a contour C in the fourth quadrant of the complex momentum plane, Eq. (I12p can be written as 
an equation for the momentum-space wave function ty(k) — 



dk'k 



>1J2 



2^ 4~ Km (kg ; A na/ 



\k')^{k')^E na ^{k). 



(14) 



In this paper the contour C is chosen to be made out of three straight-line segments Ci.2,3, C = 
C\ + C-2 + C3. Segment C\ extends from ko — to k\ = k\ r + iku, segment C2 from k\ to k2 = fc2r, and 
segment C3 from ki to k% — k max , where k max > ki r > k\ r > > kn > —k\ r are real numbers. Since 
the sp set must be finite, the contour integral along C is performed up to a cutoff k max and discretized 
with a quadrature method. In this case, k max must be chosen large enough such that all low-energy 
physics below A na is taken into account. Here, we typically chose k max ~ 3/l„ Q . Had we chosen a 
sharp regulator for F(x), k max would have been such that k max = A na since in that case V na (ko, A na ) 
would have vanished for momentum above A na . In practice, the contour C is discretized using a Gauss- 
Legendre quadrature using iVj points for the segment Cj, for a total number N s h = iVi + N% + N3 of 
discretization points. 

If ku = 0, the contour is along the real axis, and the solutions of Eq. (fl4l consist of bound states 
and scattering states. If ku ^ 0, solutions consist instead of bound states, resonant states located 
above the contour, and complex-scattering scattering states along the contour (l7l |. In order to include 
the P3/2 resonance, we take k\ r = 0.18 fm -1 , ku = —0.08 fm -1 , and k,2 r — 0.5 fm -1 . We show in 
Fig. [T]the position of the p 3 / 2 resonance as a function of A na . For each value of A na the Schrodinger 
equation flT4"|) is solved with the LO potential along the complex contour C described above. Results 
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Fig. 1 Position of the lowest p 3 / 2 not resonance as a function of the cutoff A na . Points at the top (bottom) 
are the results for the real (imaginary) momentum of the resonance coming from the numerical solution of the 
Schrodinger equation described in the text. The dashed lines are the corresponding momenta obtained from 
the effective range expansion with empirical values for the scattering volume and effective momentum. 



for the energy of the resonance are independent of the choice of the contour as long as it goes below 
the resonance and as long as the discretization is precise enough. For instance, for A na — 3.1 fm -1 we 
use k max — 10 fm -1 , N\ = 10, N 2 = 10, and -/V 3 = 15. For comparison, k res obtained directly from the 
ERE is also shown. As it can been seen from the figure, as the cutoff A na increases the position of the 
resonance quickly converges to k res . 

Since the 6 He ground state is bound, it is equivalent to use for the sp states either a set of shells 
located along the real continuum axis, or a complex-continuum set of states along a complex contour 
C along with the resonance. In the three-body calculations presented below we have used shells 
located along the real energy axis, that is, we have taken ku = 0. For instance, for A na — 7.1 fin , 
we use ki r = 5.0 fm -1 , k 2r = 12.0 fm -1 , and k rnax = 21.0 fm -1 , with Ni = N 2 = N 3 = 30. We could 
certainly decrease the number of shells to reach the same precision in the value of the 6 He ground 
state. We have not studied in detail what would be the smallest admissible number of points. Indeed, 
for a system made of three particles this is not of a great importance, since the diagonalization of the 
non-symmetric Hamiltonian matrix can be performed rather fast. 

From the set of eigenstates of the Schrodinger equation (|12p a sp basis is generated. For an energy- 
independent potential, a resonant state \& r es) above, and scattering states along, C satisfy the usual 
Berggren relation [I?} , 

\<Pres) (*res\ + j c dkk 2 \<P(k))(<f(k)\ = 1, (15) 

where the bra (<^| conjugate to the ket |<P) is such that (&\r) — (r\^). However, the potential 
V na (ko; A na ) being energy-dependent, the eigenstates of Eq. ([13)1 are not orthogonal and Eq. (fT5j) 
does not hold. One then has to consider an extra step to generate a basis. This is achieved, after 
having discretized the contour C, by solving the matrix equation 

f> t )<^ oc | = l, (16) 
i=l 

where is one of the N s h discrete sp eigenstates of the potential, Eq. (fl2l) . and (fj | its bi-orthogonal 
complement. By construction, 

(^ oc |^i>=%. (17) 

A complication is that for sufficiently large values of the cutoff, A na > A^ ~ 1.8 fm -1 , the potential 
V na (k ; A na ) supports a bound state \&b)- At A na — 1.8 fm -1 the energy of this bound state is 
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Eb = —20.941 MeV, that is, outside the range of validity of our EFT approach. As a consequence, 
we do not want to include it in the valence space. From the practical point of view, we first tried to 
construct the bi-orthogonal basis by including the bound state in the bi-orthogonalization procedure, 
and then omitting it when constructing the many-body basis to solve the three-body problem. This 
procedure turned out to give rather peculiar results in the sense that the value obtained for the 6 He 
ground state displayed a discontinuous behavior as A na varied from values below to above At,. We 
suspect that despite the fact that the bound state is not included in the basis, it has an indirect effect, 
for it is present during the phase of construction of the sp basis according to Eq. (TTBl) . One has then 
to figure out another way to generate the sp shells when a deep bound state is present. 

For A na > ylb, we generate the sp basis by converting the energy dependence of V na (ko; A na ) into 
momentum dependence by introducing an energy- independent potential V 7 ' la (A na ) that reproduces the 
half-on-shell T matrix [23j, 

(k'\V^ a (A na )\V) = (k'\V na (k ;A na )\V), (18) 

where ko is obtained from the energy of the H sp eigenstate For each discretized value k[ along the 
contour C, we solve Eq. (|18|) in order to generate V^ Q (fc^, kj\ A na ) without considering the bound state 
Itf'b). For each k[ we have N s h unknowns, V^ a (/^, kf, A na ) with j = 1, . . . , N s h, and N s h — 1 equations, 

(kl\VUAn a m) = (k'AV^k^KoWj), j = .,N ah , J ? b. (19) 
In order to solve this linear system we impose the condition 

(k'\VUAna)\k Nsh )=0, (20) 

with |fcjv s)1 ) being the state with the largest momentum on the contour C. This leads to a small error, 
since at such high momentum, fc/v sfe ~ k max ~ %A not , the influence of the potential is negligible. 
Moreover this error can be made arbitrarily small by increasing k max . The potential V^ a {A na ) is non- 
Hermitian and has right eigenvectors and left eigenvectors (^ e ^'|. The right eigenvectors are by 
construction the eigenvectors of the original energy-dependent potential V na (ko;A na ), and we now 
have the following completeness relation: 

Y^m& eft \ = i. (21) 

For A n - a < A[>, the two previous procedures for constructing the sp basis are completely identical, 
the left eigenvectors ("^ | of V^ ia (A na ) obtained with the second method being equal to the bi- 
orthogonal complement states {fii° c \ obtained in the first method by solving Eq. (|T5|) . 

From the sp basis, we construct the antisymmetrized three-body basis states coupled to good 
total angular momentum J, \{^i, 1 ^j)f < •), which are eigenstates of the Hamiltonian Hl P + Hif with 
eigenvalues Ei + Ej : 

(H? + H?) m, ^)f<;> = + E 3 ) m,*^). (22) 

The corresponding bi-orthogonal complement is ((l^ e ,<^- Yi<j\- 

The interaction V nn {A nn ) is defined in terms of relative coordinates between the two neutrons. Since 
our Hamiltonian was written in terms of na coordinates, a transformation is necessary to express the 
matrix elements V nn (A nn ) in the shell model basis #j )!<,•)■ For this purpose we use an expansion 
on a set of harmonic-oscillator (HO) wave functions, as in Ref. 16]. That is, we project V nn (A nn ) on a 
HO set \ab), where a and b label sp states of HOs in the na coordinate, and consider the nn interaction 

C(4n) = H \ab)(ab\V nn (A nn )\cd)(cd\, (23) 

a<b c<d 

where the restriction in the sum is due to the antisymmetry of the two-neutron state. Using Moshinsky 
transformations [24j, one can easily calculate (ab\V\cd). Results for the three-body energy are inde- 
pendent of the values of the HO frequency, as long as enough HO states are included in the expansion. 
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Fig. 2 Ground-state energy of 6 He from the LO two-body potentials V„ a (ko; A na ) and V n „(A nn ). For each 
value of A nn the cutoff A na is varied. 



4 Results for the 6 He Ground State 

The ground state of 6 He is coupled to J n = + and the three-body basis states are constructed from 
the sp states of the V na (ko; A na ) potential as described in the previous section. At LO, only p 3 / 2 shells 
are included in the valence space and, as a consequence, all matrix elements of the recoil term in Eq. 
(|TT|) vanish. For each value of A na the coupling constants A(A na ) and B(A na ) are fixed such that 
the ERE in the p 3 « channel truncated at the level of the effective "range" is reproduced. Similarly, 
Co(A nn ) is fixed such that the nn scattering length is reproduced. 

Figure [3 shows the energy E nna of the ground state in 6 He for different values of A na and A nn . For 
each value of A nn the cutoff A na is increased. One can see that the energy initially quickly decreases, 
then slowly rises. For A nn = 1.6 fm -1 , for example, E nna goes from —0.034 MeV for A na = 2.1 fm -1 
to -0.475 MeV for A na = 6.1 fin -1 , then to -0.400 MeV for A na = 12.1 fm -1 . As A nn increases, the 
initial decrease becomes steeper, and the increase is postponed to higher values of A na . For instance, 
at An„ = 2.5 fin -1 , the energy goes from -0.182 MeV to -2.251 MeV to -2.524 MeV in the same 
range of A na values. 

This behavior can be understood from the qualitative renormalization features of the system. As 
A na increases, the phase space of the three-body system increases, the attractive nn interaction is 
better resolved, and the binding energy increases. This is consistent with the pattern observed in Ref. 
[25l | for the energy of a three-fermion system interacting via a two-body force constructed with EFT 
at LO. In that case, for a fixed cutoff of the two-body interaction, the total energy of the system 
decreases as the size of the model space increased. As A nn increases, presumably more correlations are 
cut off for too small a value of A nai generating the faster decrease. However, there is also a residual 
dependence on A na from V na (ko; A na ). Even though the potential has been properly renormalized, 
that is, the coupling constants A(A na ) and B(A na ) have been fixed so that the truncated ERE is 
reproduced, there still is a dependence for finite values of the cutoff, as seen in Fig. [TJ The energy of 
the p 3/2 resonance goes from k = 0.7714 - 0.2947i MeV to k = 0.7696 - 0.2896i MeV when A na goes 
from 6.1 fm -1 to 12.1 fm -1 . This means that, as A na is varied within this range, there is a variation 
~ 0.005 MeV, or about 7%, in the norm of the energies of the p 3 / 2 resonance, which is consistent with 
a variation of about 15% in the three-body energy in the same range — for example a variation of 
~ 0.075 MeV for A nn = 1.6 fm" 1 . 

One can clearly see from Fig. [2] that as the cutoffs A nn and A na are increased, the energy decreases 
without reaching a stabilized value. To stress this fact, in Fig. [3] we plot the 6 He ground-state energy 
as function of A nn = A na . We have checked that the results are similar if other relations are assumed 
between A nn and A na , for example, if we take the minimum energy for each A nni which is equivalent 



Fig. 3 Ground-state energy of 6 He from the LO two-body potentials V na (ko; A na ) and Vn n (A nn ), for A nn — 
A na . 

to choosing A na large enough so that all correlations of the nn interaction have been resolved by the 
three-body system. 

The nearly linear dive of the ground state seen in Fig. [3] is reminiscent of the behavior observed 
with LO two-body forces in EFTs for systems of three bosons or three-or-more-component fermions 
[ldj . There, the dive is even faster, more like quadratic in the cutoff, stemming from the strong s-wave 
interactions among the three pairs. In either case, what we see is a collapse of the ground state under 
short-range two-body interactions similar to the one first observed by Thomas [9J. It is an indication 
that the three-body problem has not been properly renormalized with only two-body interactions [To| . 
The cutoff dependence is not decreasing as the cutoffs increase, as one would expected from residual 
cutoff dependence in a renormalized system that has been truncated correctly, but is instead increasing 
with positive powers of the cutoffs. 

The solution to this problem has to be found outside the two-body subsystems, which are perfectly 
well defined and well described by the EFT. We thus add a three-body force to renormalize the three- 
body problem. An s-wave three-body force does not have any impact at LO on the structure of 6 He 
since the na subsystems are in a relative p 3 / 2 wave. The lowest-derivative three-body force that does 
not vanish in the channel of interest can be written, in the coordinates we are using, as 

V nna (k[, k' 2 , k u k 2 ) = D(A nna ) k'^k'^ F {kf/A 2 nna ) F (k'i/Al na ) F (k\j A 2 nna ) F (k\j ' A 2 nno ) , (24) 

with k'i (kj) the outgoing (incoming) momentum for the ith na subsystem. Here A nna is a three-body 
cutoff and D is a low-energy coupling constant with dimensions of mass -9 , whose dependence on A nna 
is adjusted so that three-body observables be (nearly) cutoff independent. 

Here for simplicity we take A nna = A nn = A na . We find that we can then keep the 6 He ground-state 
energy E nna constant as the cutoff is varied. We show in Fig. 0]the resulting running of the coupling 
constant D(A nna ) when the 6 He ground-state energy is fixed to its experimental value Eq 9S = —0.98 
MeV 26]. At low cutoffs, D is negative. From Fig. [3]we see that at a cutoff Aq ~ 2.9 fm -1 the energy 
calculated with only two-body forces agrees with the experimental value, so D(Ao) — 0. Above Aq, 
DiAnn^A^^ is positive and approximately constant in the region of cutoffs we could probe. We 
cannot, however, exclude a limit-cycle-like behavior at higher cutoffs, as observed in Ref. [101 ] . 

Again like for three bosons or three-or-more-component fermions [Iol |, RG invariance requires the 
three-body force to appear at LO. Naturalness together with naive dimensional analysis suggests that 
D would scale as A/~ 9 , with M a large mass scale such as the alpha-particle binding momentum or 
the pion mass. If that were the case after renormalization, the three-body force (|24[) would be a very 
high-order effect. Instead here, as for the three-nucleon system [l(|, a certain amount of fine-tuning 
is present: the low-energy scale responsible for the existence of the shallow two-nucleon 1 Sq virtual 
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Fig. 4 Dimensionless three-body coupling constant /x 4 E^ na A^ na D as function of A nna / ^ jj,E nna , when the 
6 He ground-state energy E nna is fixed at its experimental value, for A nna — A nn = A na . 



bound state and the shallow J He p 3 / 2 resonance must appear in the renormalized three-body force 
as well. The infrared enhancement of the LO two-body interactions dominates the running of the LO 
three-body force, making its effects much larger than the naturalness expectation. While in the pure 
s-wave case the enhancement is proportional to the square of the large scattering length [l(| , here it 
must be roughly the square of the large scattering volume. 

With the three-body force so determined, we have looked for other + bound states and found none 
within the cutoff range we investigated. This is perhaps not surprising. It has been argued that the 
Ehmov effect [8[ is present if both the scattering volume and the effective momentum in a system with 
pairwise p-wave interactions are large, although there is debate about whether this can be realized (l2| . 
Since r na is not particularly large, we would not expect here an Efimov tower of shallow three-body 
states anyway. 



5 Conclusions and Outlook 

In this paper we have described for the first time the ground state of 6 He using interactions derived 
from Halo/Cluster Effective Field Theory, where the alpha-particle core is treated as an explicit field 
0. The two-body not and nn interactions are of the contact type, with parameters determined from 
two-body scattering data. The three-body dynamics of the system was solved using the formalism of 
the Gamow Shell Model [l3j], where the set of single-particle states (resonant and continuum) is given 
by the na potential. We had to adapt the formalism to accommodate the energy dependence of the 
LO na EFT potential. This is also the first time the GSM has been applied to the solution of EFT. 

We have seen that, at leading order, two-body forces are not sufficient to properly renormalizc 
the three-body system, even though they provide a systematic expansion for two-body scattering 
Indeed, as the cutoffs are increased the energy of the three-body ground state does not stabilize and 
would collapse for an arbitrarily large cutoff. We have shown that the addition of a single three- 
body counterterm is enough to achieve renormalization-group invariance0, as for systems with s-wave 
interactions [l0|. We have obtained the RG running of the coupling constant by demanding that the 
binding energy be fixed at its experimental value. 

Our work paves the way for more comprehensive studies of halo nuclei with H/CEFT. For the future, 
we plan to carry out a more extensive investigation of 6 He, including higher-order corrections and 
calculation of other observables (such as the ground-state radius and the first excited-state energy) . At 
the cost of more computational resources, other members of the He isotope family could be investigated 



2 It is our understanding that the same conclusion was reached by Ji, Elster and Phillips [T^; [27| . 
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as well, along the lines of Refs. Q EH- More generally, we hope that the combination of EFT and 
GSM will prove to be a valuable tool in the study of other three-body resonant states, such as the 
Hoyle state in 12 C. 
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